Frustration, Area Law, and Interference in Quantum Spin Models 
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We study frustrated quantum systems from a quantum information perspective. Within this 
approach, we find that highly frustrated systems do not follow any general "area law" of block 
entanglement, while weakly frustrated ones have area laws similar to those of nonfrustrated systems 
away from criticality. To calculate the block entanglement in systems with degenerate ground states, 
typical in frustrated systems, we define a "cooling" procedure of the ground state manifold, and 
propose a frustration degree and a method to quantify constructive and destructive interference 
effects of entanglement. 



The study of entanglement in spin lattice models have 
provided a novel insight into the complexity of magnetic 
ordering and quantum phase transitions in a large vari- 
ety of models In particular, in bipartite partitions 
of certain many-body systems, the block entanglement 
(BE) evaluated in the (gapped) ground state (GS) of the 
system, is proportional to the area of the block boundary. 
More complex relations arise at criticality (non-gapped 
systems). These generic relations between entanglement 
and area are known under the name of area laws 0]. 

There are several intriguing questions concerning these 
relations Q , with respect to systems with unique ground 
states. Within this context, one may ask whether the 
scaling of entanglement can also help to characterize frus- 
trated many body systems Typically, frustrated sys- 
tems exhibit a large GS degeneracy, or quasi-degeneracy, 
and a rich phase diagram [5j . Frustration in spin systems 
arise when the GS spin configuration does not simulta- 
neously satisfy all the constraints imposed by the Hamil- 
tonian; it may be caused by the lattice geometry (as in 
Ising antiferromagnets (AF) on a triangular lattice), or 
by the presence of disorder (as in spin glasses). Recently, 
the relation of frustration to high-T c superconductivity, 
as well as the discovery of "exotic" frustrated phases has 
triggered a renewed interest in the subject [H, [|| 0] • Ul- 
tracold atomic gases offer unprecedented ways to control 
many-body systems, and provide a perfect playground 
to study disordered and frustrated systems. In partic- 
ular, this should allow for experimental study of effects 
discussed in this paper (cf. @, see for a review). 

The main thesis of this paper is that highly frustrated 
systems do not follow any general area law, while weakly 
frustrated systems have area laws similar to those of 
nonfrustrated systems (away from criticality). We intro- 
duce three tools to characterize frustrated systems and 
deal with the above thesis: (i) a paradigm of a certain 
"cooled state" for finding area laws for systems which 
have a degenerate ground state manifold, including frus- 
trated systems, (ii) a "frustration degree" for quantifying 
frustration in a system, and (iii) a method for quantify- 
ing constructive and destructive interference of entan- 



glement. Using these tools, we prove the above thesis 
for six prototype frustrated spin models, each having a 
different type of frustration. More specifically, we show 
that for the Ising model with long-range (LR) AF in- 
teractions (1), and the LR AF Heisenberg model (2), 
with frustration degrees close to unity, there is a com- 
plete departure from the usual area law for nonfrustrated 
systems. However, when we consider systems whose frus- 
tration degree is much smaller, we find areas laws similar 
to those of non-frustrated systems (away from critical- 
ity), i.e. the entropy of a block of characteristic length 
L, in a system of N particles in dimension D, scales as 
L D_1 , when L D <C N. Among this class, we have an- 
alyzed a two-dimensional (2D) Ji — J 2 — J3 Heisenberg 
model with resonating valence bond (RVB) states as GSs 
(3), the Shastry-Sutherland model (in 2D) (4) [l(|, the 
one-dimensional (ID) J\ — J 2 Heisenberg model at the 
Majumdar-Ghosh (MG) point (5) [HI], an d finally an 
Ising chain with a single disordered interaction (6). 

Cooling procedure. To treat the problem of degener- 
acy of the GS manifold and calculate bipartite entangle- 
ment on this manifold, we introduce a "cooling" method. 
Suppose that a system, made up of several subsystems 
(e.g. spins), is described by the Hamiltonian Ti. We 
choose some suitable initial product state \ty ) = Jj { 
(where i runs over all subsystems) and "cool" |^o) to 
below a desired energy level E, by projecting |\& ) onto 
the subspace that is spanned by all energy eigenstates of 
7i whose energies are below E. This quenching method 
is a caricature of evaporative cooling of trapped atomic 
gases, where atoms having energy above a certain value 
are removed from the trap. So, the resulting state is 

|$ E ) = (l/\/Z)P[<E]|* ), (1) 

where P[< E] = £ E .< E Pi[ E d\ with P 4 E i] bein g the 
projector onto the eigenstates of energy Ej, and Z = 
(*o|-P[< E]|* )- For an iV-spin state |$ E ), we char- 
acterize the entanglement between k spins and the rest 
of the system by the von Neumann entropy E^-.N-k = 
5(trfc|$E) ($e|) _ the unique asymptotic measure of en- 
tanglement for pure states in bipartite splittings [13] • We 
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study the scaling of BE Ek-.N-k, for large N, where the 
initial state |^o) is chosen to maximize -Ejv/2:iv/2 13 1 . 

Frustration degree. Before proceeding further, we de- 
fine a frustration degree (J 7 ) for quantum spin mod- 
els. Frustration is a classical concept, and thus con- 
nected to the classical configurations of the quantum 
system. Given a Hamiltonian Ti., and a corresponding 
ground state \Q), we replace the one-body, two-body, 
. . . terms in H by the corresponding Ising ones (i.e. 
of, oferj, etc.) to obtain H x , discarding any constant 
term. Here o~ a (a = x,y, z) are the Pauli operators, and 
a = (a x ,a y ,a z ). We find the terms Ti) (Ti l nf ) of Ti 1 
that gives rise to positive (negative or zero) energies in 
\G). For a given \Q), Ti 1 = £ fc U ) + E; Kf We define 
the frustration degree of the system described by H as 

T = Av „ u r \T_,i lru , where the average (Av) is taken 

over all ground states of Ti. Note that the denominator 
cannot vanish. In the case of isotropic Hciscnberg •CTj) 
and XY (of ctJ + of 0jf) interactions, to find T , one may 
replace and (of, of) by SU(2) and SO(2) symmet- 
ric classical vectors, of appropriate lengths, respectively. 
With these concepts at hand, we move to study the dif- 
ferent cases. 

Case 1: Ising model with AF LR interactions: 
Ising gas. The governing Hamiltonian, for 2m 
spins, is H^iX) = (J/2m)(S - 2m\) 2 = 
(J/2m)(2E?j =1(<>i) of | - 4mA E^i< + const.), 
where S = X^of, < A < 1. In the frustrated case 
(J > 0), T = (1 + 2A - A 2 - l/m)/(l + A) 2 . For 
large systems, with A = 0, T « 1. For the initial state 
l*o> = II>|0) + with a/5 ^ (with |0) and |1) 

being the eigenstates of a z with eigenvalues ±1, respec- 
tively), the cooled state (which is independent of a, 0) is 
|$ i ^™ 9 (A)) = superposition of all states with m(l + A) 
|0)'s and m(l — A) |l)'s, where 2mA is a positive integer 
< 2m. For arbitrary 2m and magnetization A, and for 
any bipartite splitting k : 2m — k, the reduced density 
matrix can be found analytically: for even k < m(l — A), 

fe 

Pk = y^[fc,i][2m -k,mn- i] \ Wk-i) (W k -i\ /[2ro,rop], 
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(g-6)!6! ' and I^J') d en °tes the 



where (j, = A + 1, [a, 6] 
normalized equal superposition state of j|l)'s and rest 
|0)'s, i.e. the (generalized) "W" state. Any other par- 
tition can be found similarly. For a given k, Ek-2m-k 
increases (and converges to a constant) with 2m. For 
a fixed m, using Stirling approximation, E^m-k ~ 
-Eto e * lo g 2 e*, with a = ^(1 + A)''(l - \) k - r [k,i], 
for large m and k, such that k -C m. In the same limit, 
Ek-.2m-k ~ I log 2 [(1 — A 2 )fc] . The leading term to the 
scaling is \ log 2 k, independently of A - a logarithmic di- 
vergence (LD) of BE for all frustrated AF LR Ising gases. 
Since this system is effectively in "infinite" dimensions, 
there is a clear departure from any "area law" (which 



could be k 1 ^ 1 / with D — > oo). Similar divergence is 
also obtained even if the quenching is performed over a 
few low lying energy manifolds. In contrast, for the ferro- 
magnetic (J < 0) nonfrustrated LR Ising model (say, at 
A = 0), the state after cooling is (|0® 2m ) + |l® 2m ))/V2, 
and has BE £Jfc : 2 m _fc=l, irrespective of k and 2m. 

Case 2: Heisenberg model with AF LR interactions: 
Heisenberg gas. The Hamiltonian, for 2m spins, is H^r — 
{J /2m) Y^i j=i <?i • °j" {J > 0). The ground state mani- 
fold can be defined by labeling the 2m spins as "black" 
(call them B) and the remaining as "white" (call them 
W), and considering all possible "coverings" by singlets 
(valence bonds) from black to white spins. Given that 
a spin- 1/2 particle has only two orthogonal states, an 
upper bound for BE in the cooled state, for an initial 
state with two states, is a reasonable approximate up- 
per bound. Cooling the initial state j*?) = Iliee l^ 8 ), 
rijgw |^ W )-; we obtain j^^), which is an RVB formed 
by an equal superposition of all the singlet coverings. 
Consider a "system" S of k <C 2m spins (containing, 
say, b black spins and w white ones), and the remain- 
ing "environment" £. The orthonormal states of S in a 
S : £ Schmidt decomposition must be fully symmetric 
under independent permutations in the subsets of black 
and white spins. So the basis of symmetric spin multi- 
plets \b/2,mb) BnS \w/2,m w } wnS will span them, where 
m c = — |, |, c = b, w, so that there are (b + l)(w + 1) 
symmetric basis states. The Schmidt decomposition can- 
not have more terms than the number of symmetric basis 
states, so that Ek-.im-k < l°g2(^ + l)(w + 1). No matter 
what is b or w, any area law of the form E ~ k x ~ x l D for 
any D (including D = oo) is out of question. In the case 
when w = k, b — (or w = 0, b = k), the cooled state is 

fc/2 m/2 

£ |fc/2,MW J2 (-l) (p - m/2) C^/V^TI 

M=-k/2 p=-m/2 

x|(m - k)/2,p- M) m e\m/2, -p) w , 

where (C™m) 2 = [k,M + k/2][m — k,p — M + (m — 
k)/2]/[m,p+m/2] defines the Clebsch-Gordon coefficient 
(C™£ = unless < b < a in all [a, b] involved), whence 
one can show, with some algebra, that the upper estimate 
is saturated. 

Case 3: 2D J\ — J2 — J3 Heisenberg model with RVB 
ground states. Consider a (2D) square lattice of size 
2m x 2m with periodic boundary conditions (PBC), each 
site being occupied by a spin-1/2 particle. It is governed 
by a two-body J\ — J2 — J3 AF Heisenberg Hamilto- 
nian, with exchange constants J\ for nearest neighbors 
(NN) and on diagonals in the plaquettes, J2 for inter- 
plaquette NNs and diagonals, and J3 for next-nearest 
neighbors (NNN) and knight's-move-distance-away spins 
on horizontal and vertical ladders formed by the plaque- 
ttes 14 , [ll| . A plaquette is a square of four neighboring 
spins. Realization of such a "quadrumerized" lattice is 
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being hotly pursued in current experiments with atoms in 
optical lattices (see (Hf). For certain choices of the rela- 
tive strengths of the couplings, it is reasonable to assume 
14 , |l5[ that the GS configuration is formed by either 
two horitzontal singlets \HH) or two vertical ones \VV) 
in each of the plaquettes, with a fixed density, d = s/m 2 , 
of s vertical singlets in the whole state of to 2 plaquettes. 
A superposition of such states is an example of an RVB 
state, whereas ordered configurations of dimers (singlets) 
correspond to valence bond solids (Peierls order). If d is 
not too close to or 1, the dimension of the GS subspace 
grows exponentially with system size: for d = 1/2, the de- 

2 

generacy scales as 2 m . Note that in contrast to Cases 1 
and 2, the model here has rather "medium" -range inter- 
actions, and the frustration degree will correspondingly 
be smaller. For instance, for J 2 = Jz = T = 1/2. 
Consider first the scaling E p:rest , for p plaquettes. For 

the initial state |*^) = H^i(a\HH) + f3\VV)) p (unnor- 
malized, and af3 ^ 0), the cooled state (which is inde- 
pendent of a, P) is 



min(s,fe) min(s — l,m 2 — k) 



I* 



RVB/ 



E E I'W 3 '^ 



1=0 r=0 



x \r)s \J 3 r [to 2 — k, r] [to 2 — I — r, s — I — r], 

where \l)s (\f)s) ar e orthonormal states of a "system" 
S ("environment" £ ) of k (to 2 — k) plaquettes. When 
to 2 — > oo, the upper limit min(s — /,to 2 — k) is l- 
independent, and by Stirling approximation, ln[m 2 — t, s— 
t] w (to 2 — t) ln(m 2 — t) — (s — t) ln(s — t) plus functions 
of to 2 and s, where t = I + r. Linearizing the logarithm 
around the mean t of t (in \fpgyg)), ana - exponentiating 
back, we obtain a product: [to 2 — t, s — t] ss q l q r , with 
q = (s-I)/(m 2 -i) m + + 12d(l - d))/(6(l-d)). 
The approximations are systematic, as t < \/3fc/2, and 
dispersion of t ~ Vl. So the states of £ and S can 
be well approximated as product states, which further 
implies that \^jivb)j f° r an Y can be approximated 
by product over plaquettes, where each plaquette is in 
the same pure state - a "mean-field picture" is valid, 
for large to 2 . Consequently, the entanglement of a sys- 
tem S of size k, plus a boundary B that intersects 
A plaquettes, to the rest of the lattice, will scale as 
hE p i(d) + vE p i(l — d), where the boundary intersects h 
(v) plaquettes horizontally (vertically), h + v = A, and 

E pl (d) = log 2 (1 + 3q 2 ) - log 2 q 2 . 

For the most general initial state, \^q s ) = 

IlSi(llMWhfe>llk>)p, denoting P^ p ° for the pro- 
jector onto the space spanned by {\HH) 7 \VV}} 
at the p th plaquette, and {|Gj)} for the ground 
states, the cooled state is (Y^,i\ G i)( G i\)\^o s ) = 

(EJG l )(G 4 |)n™>^ ( W s ) = (EijGWGi|)|*?>, 
so that the same area law holds. Surprisingly therefore, 
the area law depends on the path of the boundary (and 



not only on its length): it depends on h, even for fixed 
A. However, the area law is linear in this 2D system, and 
supports our thesis for weakly frustrated systems. 

Case 4 : Shastry- Sutherland model. It is a 2D NN an- 
tiferromagnet (with Heisenbcrg couplings of strength J\) 
on a square lattice with additional Heisenberg in- 

teractions of strength J 2 on the diagonals (2i, 2j) «-» (2i+ 
Mj+l) and (2i,2j+l) (2i-l, 2j+2), withPBC. The 
ground state of this model, for Ji/Jo < 1/2, is a product 
of dimers along the J 2 diagonals [lflj. For Ji/J 2 < 1/2, 
by going to the Ising limit, the degeneracy becomes ex- 
ponentially large, and T » 1/(1 + (1/2)(J 2 /Ji)) < 1/2, 
for large system size. Furthermore, the area law is linear, 
and depends on the path of the boundary. 

Case 5: ID J\-Ji Heisenberg model. Consider a ID 
system with PBC and 2to sites, having NN and NNN 
Heisenberg couplings of strengths J\ and J 2 respectively. 
For Ji = 2J 2 > (MG model [13), the GSs are 
\Gmg) = Uti (|0>2i|l) 2i ±i - |l> 2i |0>2i±i) /V2. At the 
MG point, T = 1/2. Note however that if we replace &i 
by a 3D classical vector, the frustration degree is zero. 
Using the initial state \^L^ |0)2i-i|l)2»|^)am-i|^)2mj 
where \4>) is an arbitrary (qubit) state, we obtain a lower 
bound: E^. rest > 2 or 1, for even or odd k. Any cooled 
state is of the form |$mg) = H^mg) + &|C?mg)' whence, 
after writing in Schmidt dccompostion, an upper bound 
reads: Ek :rest < log 2 5 or log 2 3, for even or odd k. Nu- 
merical simulations at the MG point, show that e.g. for 
2m = 8, E 4:rest « 2.314 (cf. 17J), and BE after cooling 
converges with k. Therefore, the area law is a constant 
for this ID system, in support of our thesis. Note that the 
above method of finding bounds can potentially be used 
in other models whose GS space is made up of dimers. 

Case 6: Ising spin chain with NN interactions. 
The Hamiltonian, of 2to spins (with PBC), is 



H 



Isijig 



NN - E(ij> J ij a i a j> with a11 \ J ij\ = J i and 

all except one are negative. T = 1/(2to — 1). 

The initial state \V$ N ) = Ui((\0) + 

(cf. [3), gives the cooled state \<&nn 9 ) = 

1 J^™ [(|0®( 2m ~ fe )l® fe ) + |l®(2m-fc)Q®fc^ _|_ 
n^®(fe+l)Q®(2m-fc-l)\ _|_ |Q(8(fc+l)jl8(2m-fc-l)\\l p or 

any fc, Ek;2m—k decreases with 2m, contrary to previous 
cases. Moreover, Ek-.2m-k, for a fixed 2to is a constant 
with k - in support of our thesis. 

Constructive and destructive interferences. The above 
studies allows us to identify an interesting interplay be- 
tween frustration and interference of entanglement. Be- 
ginning with Case 3, for a fixed density d, the lattice 
is in the state IV^kb)) which is an equal superposition 
over all states that have s vertical plaquettes and to 2 — s 
horizontal ones. If we choose one of the states in this 
superposition, calculate the entropy of the inner area 
(S U B) in the chosen state, and then average this en- 
tropy over all states in the superposition, the average 
entropy is E(d) = 2hd + 2v(l — d) . For a square B, 
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Density d Density d 

FIG. 1: Constructive and destructive interference. (Left) 
E/E for a square B, so that v — h. (Right) E/E for a horizon- 
tal (vertical) boundary, so that v = (h — 0). The vertical 
axes are of E/E. 

the interference is destructive (i.e. E/E < 1) for low 
and high densities, but otherwise constructive, and for a 
purely horizontal (or vertical) B, interference can wipe 
out the entanglement completely, at low densities (Fig. 
[T]). For a square B, E/E is maximal for d — 1/2. Cor- 
respondingly, in the AF LR Ising gas, the highest BE 
scaling is for A = 0, when there are an equal number |0)s 
and |l)s in the components of the cooled state, and T is 
also maximal exactly at A = 0. We predict that a parent 
Hamiltonian that describes the ground states in Case 3 
for fixed density, will be maximally frustrated at d = 1/2. 

For a similar study in the AF LR Ising gas (Case 1), 
we normalize the entanglement scaling in the frustrated 
system, i.e. J > 0, with the one in the nonfrustrated one 
(J < 0). Since the latter is unity, the normalized entan- 
glement, to leading order, is \ log 2 k, for large blocks of 
size k, large system size 2m (3> k), and for all A. Thus we 
obtain constructive interference of BE for all A, with re- 
spect to the nonfrustrated situation. For the MG model, 
let us once again choose any one of the dimer ground 
states, find the entanglement in a bipartite split, and av- 
erage over the two ground states. The result is unity, for 
any split, so that we again have constructive interference 
of entanglement. For the system described in Case 6, we 
have marginally constructive interference, compared to 
the nonfrustrated case (all Jy < 0). 

In Case 2, we have a rather remarkable example of de- 
structive interference, because when 2m 3> k, Ek-.rest ~ k 
for most elements in the superposition forming the cooled 
state, while the latter has only logarithmic scaling at 
most. Superposition can therefore give rise to qualita- 
tive changes in scaling. 

Summary. Firstly, our studies show that the dimension 
of the GS manifold does not provide a "good" character- 
ization of frustration. This can be seen by comparing the 
Ji — J2 model at the MG point with the ID Ising model 
with a single frustrated bond. Secondly, we found that 
trying to confer an "area law" on a frustrated system can 
have surprising consequences, such as logarithmic diver- 
gence of BE in an effectively infinite dimensional system, 
and dependence of BE on the shape (and not only the Eu- 
clidean area) of the boundary. Interestingly, the seeming 



independence of area law on dimension, in the long-range 
Ising model, gives rise to the possibility of applying den- 
sity matrix renormalization group techniques to complex 
systems with long-range interactions (cf. Our re- 

sults indicate that while weakly frustrated systems follow 
the usual area law known in the literature, strongly frus- 
trated systems will each have their own area law. Finally, 
we have introduced a cooling procedure to study BE in 
degenerate ground state manifolds, a frustration degree, 
and a method to quantify constructive and destructive 
interference of entanglement. 

We acknowledge comments of H.-U. Everts, R. Fazio, 
T. Vekua, and support of research project of Polish 
Government scientific funds (2005-08), ESF QUDEDIS, 
Marie Curie ATK project COCOS (MTKD-CT-2004- 
517186), Spanish MEC (FIS-2005-04627/01368, Con- 
solider Project QOIT, & Ramon y Cajal), & EU IP 
SCALA. 



[1] L. Amico et al , |quant-ph/0703044| 

[2] G. Vidal et al, Phys. Rev. Lett. 90, 227902 (2003); G. 
Refael and J.E. Moore, ibid. 93, 260602 (2004); V.E. Ko- 
repin, ibid. 92, 096402 (2004); K. Audenaert et al, Phys. 
Rev. A 66, 042327 (2002); B.-Q. Jin and V.E. Korepin, 
ibid. 69, 062314 (2004); B.-Q. Jin and V.E. Korepin, J. 
Stat. Phys. 116, 79 (2004); P. Calabrese and J. Cardy, 
J. Stat. Mech. P06002 (2004); I. Peschel, ibid. P12005 
(2004); A.R. Its et al, J. Phys. A 38, 2975 (2005); M. 
Cramer and J. Eisert, New J. Phys. 8, 71 (2006); A. Ki- 
taev and J. Preskill, Phys. Rev. Lett. 96, 110 404 (2006); 
A. Riera and J.I. Latorre, quant-ph/0605112 

[3] J.I. Latorre et al, Phys. Rev. A 71, 064101 (2005); 
M.M. Wolf, Phys. Rev. Lett. 96, 010404 (2006); D. 
Gio ev and I. Klich, ib id. 100503 (2006 ); T. Barthel et 
al, |cond-mat706 02077; W. Li et al, |quant-ph706 02094 
T. Barthel et al, Phys. Rev. Lett. 97, 220402 (2006); j! 
Vidal et al, J. Stat. Mech. P01015 (2007). 

[4] This question was recently adressed for resonating va- 
lence bond states in F. Alet et al, cond-mat/0703027 

[5] G. Misguich a nd C. L huillier, |cond-ma t/0310405 
C. Lhuillier, |cond-mat/0502464| F. ATet" et al, 
|cond-mat/0511516| ~ 

[6] M. Rasolt and Z. Tesanovic, Rev. Mod. Phys. 64, 709 
(1992); M. Sigrist and T.M. Rice, ibid. 67, 503 (1995). 

[7] R. Melzi et al, Ph ys. Rev. Lett. 85 , 1318 (2000). 

[8] J. Billy et al , larXiv:0804.T62ll G. Roati et al, 
larXiv:0804.2609l 

[9] M. Lewenstein et al, Adv. in Phys. 56, 243 (2007). 
[10] B.S. Shastry and B. Sutherland, Physica (Amsterdam) 

108B+C, 1069 (1981). 
[11] C.K. Majumdar and D.K. Ghosh, J. Math. Phys. 10, 
1388 (1969). 

[12] C.H. Bennett et al, Phys. Rev. A 53, 2046 (1996). 

[13] If the scaling is not unique after the maximization, we 
maximize entanglement in the partition (N/2) — 1 : 
(N/2) + 1, etc., until a unique scaling is obtained. 

[14] B. Kumar, Phys. Rev. B 66, 024406 (2002). 

[15] M. Mambrini et al, Phys. Rev. B 74, 144422 (2006). 



5 



[16] B. Paredes and I. Bloch, arXiv:0711.3796 
[17] R.W. Chhajlany et al, |cond-mat/0610379| 
[18] H.J. Briegel and R. Raussendorf, Phys. Rev. Lett. 86, 
910 (2001). 



[19] J. Rodrfguez-Laguna and G.E. Santoro, 
|cond-mat/0610661| A. Ferraro et al, |quant-ph/0701009| 



